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Abstract 

A striking correlation has recently been observed between global cloud cover and 
the flux of incident cosmic rays. The effect of natural variations in the cosmic ray 
flux is large, causing estimated changes in the Earth's energy radiation balance that 
are comparable to those attributed to greenhouse gases from the burning of fossil 
fuels since the Industrial Revolution. However a direct link between cosmic rays 
and cloud formation has not been unambiguously established. We therefore pro- 
pose to experimentally measure cloud (water droplet) formation under controlled 
conditions in a test beam at CERN with a CLOUD chamber, duplicating the condi- 
tions prevailing in the troposphere. These data, which have never been previously 
obtained, will allow a detailed understanding of the possible effects of cosmic rays 
on clouds and confirm, or otherwise, a direct link between cosmic rays, global cloud 
cover and the Earth's climate. The measurements will, in turn, allow more reliable 
calculations to be made of the residual effect on global temperatures of the burning 
of fossil fuels, an issue of profound importance to society. Furthermore, light radio- 
isotope records indicate a correlation has existed between global climate and the 
cosmic ray flux extending back over the present inter-glacial and perhaps earlier. 
This suggests it may eventually become possible to make long-term (10-1,000 year) 
predictions of changes in the Earth's climate, provided a deeper understanding can 
be achieved of the "geomagnetic climate" of the Sun and Earth that modulates the 
cosmic-ray flux. 
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1 Introduction 



Global warming is a major concern of the world, with its potentially devastating effects 
on coastal settlements and world agriculture. The steep rise in greenhouse gas emissions 
since the Industrial Revolution has increased the CO2 concentration in the atmosphere 
by 30% (Fig. |lD. This is widely believed to be the dominant cause of the observed rise of 
about 0.6°C in the global mean surface temperature during this period (Fig. |^) [|T|. 

A small systematic rise or fall in the global temperature is caused by a net imbalance 
("forcing") in the Earth's energy radiation budget (Fig. The radiative forcing caused 
by the increase in the CO2 fraction since 1750 is estimated to be 1.5 Wm~^ |l[], compared 
with the global average incoming solar radiation of 342 Wm~^, i.e. an imbalance of 
only 0.4%. After including the effects of all greenhouse gases (+2.45 Wm~^), aerosolsQ 
(-0.5 Wm~^) and their effects on clouds (-0.75 Wm~^), the present net radiative forcing 
from mankind is estimated to be 1.2 Wm~^. 

The climate models |l[] upon which the predictions of greenhouse warming depend 
have gradually improved as new effects and better data have been incorporated. They 
now provide a reasonable representation of the observed variation in global temperature 
over the last century (Fig. H). However they remain subject to significant uncertainties, 
both from unknown processes and also from the unknown effects of feedback mechanisms 
which may amplify or damp a warming trend. Important unknown effects include those 
due to increased evaporation of water (which may cool due to an increased cloudiness 
or warm due to an increased greenhouse effect), shifts in the amount of CO2 dissolved 
in the oceans, changes in the ocean currents and changes in the polar ice. Given the 
uncertainties in these models, many scientists consider there is no convincing evidence at 
present to demonstrate that a greenhouse warming of the planet is in progress, let alone 
that the models can reliably predict future temperature changes (the predicted rise by 
the end of the next century is 1.3-2.5°C; Fig. ^). 

Nonetheless it is in this climate of scientific uncertainty that major political decisions 
on greenhouse gas emissions are presently being made (Earth Summit in Rio de Janeiro, 
1992, and UN Climate Convention in Kyoto, 1997) that will have a profound effect on 
the economic development of both the developed and the developing countries. The need 
for such political decisions to be based on sound scientific grounds is self evident, and a 
major world-wide research effort on climate change is underway. 

An unexpected and important factor in climate change has recently been discovered: 
cloud cover seems to follow natural variations in the incident cosmic ray flux p, ^ . How- 
ever a direct link between cosmic rays and cloud formation has not been unambiguously 
established and, moreover, the microphysical mechanism is poorly understood. This 
presents particle physics with a unique opportunity to make a major contribution to 
the problem of global warming by confirming, or otherwise, this link under controlled 
conditions in a test beam. The purpose of this paper is to provide an initial outline of the 
motivation and conceptual design of the proposed CLOUD (Cosmics Leaving Outdoor 
Droplets) detector. 

^Aerosols are 0.001-1 /zm diameter particles of liquid or solid in suspension. Atmospheric aerosols 
include dust, sea salt, soot (elemental carbon), organic compounds from biomass burning, sulphates 
(especially (NH4)2S04) from SO2, and nitrates from NO and NO2. Aerosol concentrations vary from a 
few per cm'^ in clean maritime air to 10^ cm~^ in highly-polluted city air. 
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Figure 1: Carbon dioxide concentration in the atmosphere over the past 1000 years 
The data are based on Antarctica ice core records and, since 1958, on direct measurements 
from Mauna Loa, Hawaii. Also shown are the estimated yearly emissions of CO2 from 
burning fossil fuels since the Industrial Revolution. 

2 The origins of climate change 

The Earth's climate has not always been the same; significant changes have occurred dur- 
ing recorded history and even greater changes have taken place since the human ancestral 
line began about 5 million years ago. The urgent questions facing science are why these 
natural changes took place and whether there is currently an influence of mankind on 
the climate. Only when the natural causes are understood and properly accounted for in 
climate models can we expect meaningful results on the anthropogenic influences. 

Several natural effects are well established. These include the Milankovitch cycles — 
variations in the Earth's orbit, tilt and spin with respect to the Sun — which seem to 
control the timing of glacials and interglacials. These cold and warm phases correspond 
to cycles of about 100,000 years and 10,000 years, respectively. (We are presently ap- 
proaching the end of an interglacial; sometime in the next 1000 years or so the climate 
will probably become much colder.) Other important natural sources of climate change 
are the anomalous warm Pacific current known as El Nino, and volcanoes that pour ash 
into the stratosphere, reducing the radiation from the Sun reaching the Earth's surface. A 
recent volcanic example is Mt. Pinatubo in the Philippines, which erupted in June 1991 
and caused a global cooling, by up to 0.4°C, over a period of about 3 years. In contrast, 
a natural effect that has until recently been harder to understand is the apparent link 
between solar activity — the sunspotQ cycle — and the weather. 

^Sunspots are areas of the Sun's photosphere where strong local magnetic fields (typically 2500 gauss, 
to be compared with the Earth's field of about 0.3 gauss) emerge vertically They appear dark because 
their temperature is about half of the surrounding photosphere (3,000 K compared with 5,800 K). They 
are thought to be generated by the different rotation rates of regions of the Sun: one revolution takes 
25 days at the equator and 28 days at mid-latitudes. This modifies the normal convective motions of the 
plasma and creates "knots" of strong magnetic field which penetrate the photosphere and provide escape 
paths for the plasma, forming sunspots. 
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Figure 2: Observed changes in global mean surface temperature from 1861 to 1994, com- 
pared with the climate model of the Intergovernmental Panel on Climate Control (IPCC) 
on the effects of greenhouse gas emissions |1[]. The model was run with a) greenhouse 
gases alone, b) greenhouse gases and aerosols, and c) greenhouse gases, aerosols and an 
estimate of changes in solar irradiance. The broken curves are the model predictions un- 
der several different assumptions of the climate sensitivity, AT2X , which is the change in 
global mean temperature for a doubling of the atmospheric CO2. 
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Figure 3: The Earth's energy radiation budget. The width of each band is proportional to 
the associated energy flux, which is computed as a global annual mean. Heat absorbed by 
the atmosphere produces a greenhouse effect since the radiation lost to space comes from 
the cold tops of clouds and from parts of the atmosphere much colder than the surface. 
The figure is adapted from reference |T| . 




Figure 4: Projected global mean surface temperatures from 1990 to 2100 based on the 
current IPCC climate model (Fig. |^c) for the full set of IS92 CO2 emission scenarios 
A climate sensitivity, AT2X = 2.5°C, is assumed. 
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Figure 5: Variation of the (smoothed) sunspot number since systematic (telescope) obser- 
vations began. The data show a) the period 1610-1740 0] and b) the period 1700-1995 
[the data are from NOAA — the National Oceanic and Atmospheric Administration]. 

2.1 Activity of the Sun 

The observation that warm weather seems to coincide with high sunspot counts and cool 
weather with low sunspot counts was made as long ago as two hundred years by the 
astronomer William Herschel who noticed that the price of wheat in England was lower 
when there were many sunspots, and higher when there were few. The most well-known 
example of this effect is known as the Maunder Minimum , the Little Ice Age between 
1645 and 1715 — which ironically almost exactly coincides with the reign of Louis XIV, 
le Roi Soleil, 1643-1715 — during which time there was an almost complete absence of 
sunspots (Fig. |a). During this period the River Thames in London regularly froze 
across and fairs complete with swings, sideshows and food stalls were a standard winter 
feature. 

Quantitative proof of a correlation between the Sun's activity and the Earth's temper- 
ature was presented by Friis-Christensen and Lassen |^ in 1991. They used the sunspot 




Figure 6: Variation during the period 1861-1989 of the sunspot cycle length (solid curve) 
and the temperature anomaly of the Northern Hemisphere (dashed curve) . The tem- 
perature data from the IPCC IQ. 

cycle length as a measure of the Sun's activity. The cycle length averages 11 years but 
has varied from 7 to 17 years (Fig. |^), with shorter cycle lengths corresponding to a more 
magnetically-active Sun. A remarkably close agreement was found between the sunspot 
cycle length and the change in land temperature of the Northern Hemisphere in the period 
between 1861 and 1989 (Fig. ^). The land temperature of the Northern Hemisphere was 
used to avoid the lag by several years of air temperatures over the oceans, due to their 
large heat capacity. This figure covers the period during which greenhouse gas emissions 
are presumed to have caused a global warming of about 0.6°C. Two features are of par- 
ticular note: firstly the dip between 1945 and 1970, which cannot be explained by the 
steadily rising greenhouse gas emissions but seems well-matched to a decrease in the Sun's 
activity, and secondly the close correspondence between the two curves over this entire 
period, which would seem to leave little room for an additional greenhouse gas effect. 

In the absence of sufficiently sensitive measurements, people suspected that the Sun's 
irradiance may be fluctuating over the solar cycle. However, from our knowledge of how 
the Sun produces energy, this would seem difficult. Radiation from the Sun was created in 
the core and, in a random walk with a mean free path of about 10 cm, has taken about fifty 
million years to reach the surface. This would smooth out any substantial fluctuations 
in radiated energy on a timescale significantly less than about a million years. Indeed 
the steadiness of the Sun's irradiance over a complete sunspot cycle has recently been 
confirmed by satellite measurements (Fig. 0|). The solar fiux is slightly higher at 
sunspot maximum; although sunspots are cooler and have reduced emission, this is more 
than compensated by an associated increase in bright areas known as plages and faculae. 
The mean irradiance changes by about 0.1 % from sunspot maximum to minimum which, 
if representative over a longer time interval, is too small (0.3 Wm~^, globally-averaged) 
to account for the observed changes in the Earth's temperature. 
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Figure 7: Satellite measurements of the variation of the Sun's irradiance over a 1.5 sunspot 
cycles, showing daily mean values and uncertainties . Approximately, sunspot maxima 
occurred during 1980 and 1990, and sunspot minima occurred during 1985-86 and 1996. 
The rapid fluctuations at the maxima are due to sunspots rotating into the field of view. 
The data are from the Active Cavity Irradiance Monitor (ACRIM I and II) and from the 
Nimbus 7 Earth Radiation Budget (ERB) and Earth Radiation Budget Satellite (ERBS) 
experiments. 



2.2 The missing link: cosmic rays? 

Experimental observations. It has been recognised for several decades that the cos- 
mic ray flux reaching the Earth is strongly modulated by the solar wincl0 which, in turn, 
is strongly influenced by the sunspot cycle (Fig. H). At times of high sunspot activity, the 
solar wind is stronger and this partially shields the lower-energy galactic cosmic radiation 
from penetrating the inner solar system and reaching the Earth. In addition, the lower- 
energy cosmic rays are affected by the geomagnetic field, which they must penetrate to 
reach the top of the atmosphere. The minimum vertical momentum of primary charged 
particles at the geomagnetic equator is about 15 GeV/c, decreasing to below 0.1 GeV/c 
at the geomagnetic poles. The modulation of the cosmic ray flux incident on the Earth's 
atmosphere at the higher geomagnetic latitudes is as large as 50% between sunspot max- 
ima and minima [§] . It represents one of the largest measurable effects of sunspot activity 
near the Earth's surface. 

But how could cosmic rays affect the Earth's weather? The breakthrough was made by 

■^The solar wind is a continuous outward flow of charged particles (mainly protons and electrons, with 
5% helium nuclei) from the plasma of the Sun's corona. The main sources are near the Sun's poles; other 
important sources include coronal holes and sunspots. The solar wind creates the huge heliosphere of 
the Sun that extends out to Neptune and beyond. At the Earth's orbit the solar wind has a velocity 
of 300-800 km s'^ {(3 = 0.001-0.003) and an intensity of (0.5-5)xl0^ particles cm^^ g-^ generating a 
magnetic field of about 5 x 10~^ Gauss. 
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Figure 8: Variation with time of sunspot number and cosmic ray flux, as measured by 
ground-based neutron counters. The neutron data are from the Univ. Chicago Neutron 
Monitor Stations at Chmax, Colorado (3400 m elevation; 3 GeV primary charged particle 
cutoff), Huancayo, Peru (3400 m; 13 GeV cutoff) and Haleakala, Hawaii (3030 m; 13 GeV 
cutoff). The stronger modulation of the cosmic ray flux at higher latitudes (Climax) is 
due to the lower primary cutoff energy. 
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Figure 9: Absolute percentage variation of global cloud cover (data points; left-hand scale) 
and relative percentage variation of cosmic ray flux (solid curve, normalised to May 1965; 
right-hand scale) |^. Also shown is the solar 10.7 cm microwave flux (dashed curve, 
in units of 10^^^ Wm~^Hz~^; right-hand scale). The cloud data are restricted to oceans; 
Nimbus 7 (triangles) and DMSP (diamonds) data are for the total Southern Hemisphere 
over oceans, and ISCCP (squares) data are for oceans with the tropics excluded. The 
cosmic ray data are neutron measurements from Climax (see Fig. H). All data are 
smoothed using a 12-month running mean. 
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Svensmark and Friis-Christensen in 1997 [0, |^ who showed a striking correlation between 
global cloud cover and the incident cosmic ray flux (Fig. Moreover, the cloud cover 
was found to have a poorer correlation with variations in the solar microwave flux (Fig. 
The latter is known to follow closely variations in the total solar irradiance, soft X-rays 
and ultraviolet rays. The modulation in cloud cover was found to be more pronounced at 
higher geomagnetic latitudes, consistent with the reduced shielding effect of the Earth's 
magnetic fleld at higher latitudes (see Fig. |). 

Over a sunspot cycle, the absolute variation in global cloud cover is about 3% (equiv- 
alent to 4.8% relative fraction), and the variation in the neutron flux is about 15-20% 
(for a primary charged particle cutoff of 3 GeV). The neutrons are mostly produced by 
primary hadronic interactions in the flrst 1-2 Xint of the atmosphere and therefore mea- 
sure the changes in cosmic ray intensity at altitudes above about 13 km. The primary 
cosmic radiation is about 80% protons, 15% He nuclei and 5% heavier nuclei. At sea level 
the most numerous charged particles are muons and their fluctuation is less pronounced, 
about 3% over a solar cycle since they are produced from a stiffer primary spectrum 
of cosmic radiation, which is therefore less affected by the solar wind. 



Fraction of clouds affected by cosmic rays. These data suggest that cosmic ra- 
diation is, under suitable atmospheric conditions, either extending cloud lifetimes (for 
example, by increasing the droplet number density) or inducing some additional cloud 
formation. The currently-accepted physical processes that affect the production of cloud 
condensation nuclei [1^, including both aerosols and ice particles, do not consider ion- 
induced effects. The new satellite observations suggest that the ions or radicals produced 
in the atmosphere by cosmic rays may somehow be affecting the nucleation, growth or 
activation of atmospheric aerosols, or the creation of ice particles. 

It is interesting to estimate the total cloud fraction that may be caused by cosmic 
rays. It is not yet known at which altitude the cosmic radiation may be affecting clouds, 
i.e. whether the clouds more closely follow the amplitude modulation of the neutron flux 
(20%) or of the muon flux (3%). However, making the most conservative assumption, a 
20% change in cosmic ray flux causes a 4.8% relative change in cloud cover. If the cloud- 
forming process is directly proportional to cosmic flux then this suggests about 25% of the 
Earth's clouds could be affected by cosmic rays. However, if the process is proportional 
to ion density then the fraction is twice as large, since a 20% increase in cosmic ray flux 
leads to a 10% increase in equilibrium ion density (see Section 3^). In summary, a large 
fraction of the Earth's clouds could be affected by cosmic radiation. 



Effect on the Earth's radiation budget. Thin clouds at high and middle altitudes 
cause a net warming due to an increased trapping of outgoing longwave radiation whereas 
optically thick clouds produce a net cooling due to a dominant increase in the albedo (re- 
flection) of the incoming shortwave solar radiation. At present there is no indication which 
clouds may be affected by cosmic rays, so we will assume for the present estimate that 
all are affected equally. Present estimates from the Earth Radiation Budget Experiment 
(ERBE) indicate, overall, that clouds reflect more energy than they trap, leading to a net 
cooling of 28 Wm~^ from the mean global cloud cover of 63% (Table |1] Ull]]). The 4.8% 



relative variation in cloud cover due to variations of the cosmic ray flux over a solar cycle 
therefore implies a change in the Earth's radiation budget of 1.3 Wm~^ (0.4% of the total 



incident solar radiation). This is a significant effect: in the 3-year period 1987-1990, the 
effect of the reduction in cosmic ray flux was comparable to the total estimated radiative 
forcing of 1.5 Wm^^ from the increase in CO2 concentration since 1750. 

Table 1: Global annual mean forcing due to various types of clouds, from the Earth 
Radiation Budget Experiment (ERBE) |jll|. The sign is defined so that positive forcing 
increases the net radiation budget of the Earth and leads to a warming; negative forcing 
decreases the net radiation and causes a cooling. 



Parameter 


High clouds 


Middle clouds 


Low clouds 


Total 


Thin Thick 


Thin Thick 


All 


Global fraction (%) 


10.1 8.6 


10.7 7.3 


26.6 


63.3 


Forcing (relative to clear sky): 
Albedo (SW radiation) (Wm'^) 
Outgoing LW radiation (Wm~^) 


-4.1 -15.6 
6.5 8.6 


-3.7 -9.9 
4.8 2.4 


-20.2 
3.5 


-53.5 
25.8 


Net forcing (Wm^^) 


2.4 -7.0 


1.1 -7.5 


-16.7 


-27.7 



Open questions. Despite the observation of Svensmark and Friis-Christensen (Fig. H), 
a direct link between cosmic rays and cloud formation has not been unambiguously es- 
tablished. The satellite data are a composite of several datasets which are subject to 
inter-calibration uncertainties. Furthermore the data are spatially and temporally selec- 
tive - they are restricted to oceans and exclude the tropics and polar regions, and they 
are largely daytime only. 

It has been argued that the single sunspot cycle for which there are reliable satellite 
measurements is insufficient to demonstrate a clear connection and that several more 
cycles (each of course requiring 11 years) are needed. Even if this extended study were 
to confirm the correlation between the sunspot cycle and cloud cover, it may not prove 
the missing link is cosmic rays. While cosmic rays certainly vary with the sunspot cycle, 
other factors that change during the cycle may be affecting clouds, quite independently 
of cosmic rays. For example, the cause may be variations of the solar ultra-violet (UV) 
radiation. Although the total irradiance of the Sun varies by only 0.1 % over the solar 
cycle, the variation is more pronounced in the UV. Wavelengths below 320 nm, which 
account for only about 2% of the total solar irradiance, vary by about 5% (200-320 
nm) to 50% (100-150 nm), or even more at shorter wavelengths |TB|. The UV radiation 
is absorbed by ozone high in the stratosphere which warms as a result. This has the 
potential to influence large scale tropospheric dynamics and hence, perhaps, cloudiness. 

Finally, the evidence of increased variations in cloud production at higher latitudes, 
while consistent with the cosmic ray explanation, could simply be due to the variations 
with latitude of temperature, air circulation, aerosols and other atmospheric conditions 
that affect cloud formation. 
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2.3 History of cosmic rays and climate change 

If the periodic 11-year cycles of the sunspots and the associated cosmic ray flux were the 
end of the story, then it would be of limited concern since there would be no resultant 
long-term change in the Earth's weather but simply another cyclic "seasonal" change 
(albeit with a period of 11 years). However there is clear evidence of longer-term and 
unexplained changes both in the Sun and in the Earth's magnetosphere and these, in 
turn, seem to have had long-term effects on the Earth's climate. 

2.3.1 Global warming during the present century 

Over the course of this century the cosmic ray flux has been steadily decreasing; it is 
weaker today at its maximum during the sunspot cycle than it was at its minimum 
around 1900. This has been revealed by analysis of the ^°Be concentration in a Greenland 
ice core (Fig. ^ [O])- Beryllium-10 is produced in the atmosphere by cosmic rays and, 
despite a lower production rate than ^"^C (Section |2.3.2| ) by about a factor 100 (3.5 x 10~^ 
atoms cm~^s~^), it has the advantages of settling out relatively rapidly (~ 2 years) and 
of a long half-life, 1.51 x 10^ years. A systematic decrease in the cosmic ray flux of the 
magnitude indicated by the ^°Be data would have caused a reduction in cloud cover and 
consequent warming of the Earth, sufficient to account for a large part of the observed 
rise of 0.6°C in global temperatures over this century which is presently attributed to 
anthropogenic greenhouse gases. 

The cause of this systematic decrease of the cosmic flux has been a strengthening 
of the solar wind over this century. This is revealed in measurements of the so-called 
aa geomagnetic index of short-term (up to 3-hour) variations in the geomagnetic fleld 
at the Earth's surface [0, 0, which is affected by the interactions of the solar wind 
with the Earth's magnetosphere (Fig. \n\ ). With this hindsight, the previous observation 
by Friis-Christensen and Lassen that the sunspot cycle length is closely correlated with 
global temperatures (Fig. ^ can now perhaps be understood: the sunspot cycle length is 
a direct measure of the strength of the solar wind and a better measure than the sunspot 
number itself. Shorter sunspot cycles correspond to a stronger solar wind and therefore 
to decreased cosmic flux, decreased cloud cover and increased temperatures. 

Although a comparison of global temperatures with direct measurements of cosmic 
rays can only be made since 1937, when systematic measurements of charged cosmic rays 
began with electronic detectors, the two are indeed found to correlate well (Fig. |1^ [^). 

2.3.2 Climatic change over the past millennium 

The observations of a correlation between the Sun's activity, cosmic rays and the Earth's 
climate may be extended to earlier times with ^^C data. This isotope is continuously 
formed in the atmosphere by low energy neutrons from cosmic radiation at a rate of 
about 2.5 atoms cm~^s~^, in the reaction: 

+ n + p 

The ^^C is rapidly oxidised to form ^^C02 and then decays by (3~ emission with a half life 
of (5730 ± 40) years. The turnover time of CO2 in the atmosphere is quite short, about 
4 years, mostly by absorption in the oceans and assimilation in living plants. However, 
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Concentration of ^°Be in a 300 m ice core from Greenland over the last 
The data are smoothed by an approximately 10 year running mean and 
have been shifted earlier by 2 years to account for settling time, b) The sunspot cycle 
over the same period, which shows a negative correlation with the short-term (~ 11 year) 
modulation of the ^°Be concentration. 
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Figure 11: Variation of the aa geomagnetic index over the last 130 years: a) yearly 
measurements and b) data smoothed by an 11-year running mean. The sunspot cycle 
over the same period is indicated in the upper figure and shows a positive correlation 
with the short-term (~ 11 year) modulation of the aa geomagnetic index. 
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Figure 12: Northern hemisphere mean temperature (dotted hne; left-hand scale) and 
charged cosmic ray flux (thick solid line; right-hand scale) since 1940 f^. Both curves are 
smoothed by an 11-year average. The diamonds show the annual temperatures, before 
smoothing. 



recirculation from the oceans has a smoothing effect such that changes in the ^"^C fraction 
on timescales less than a few decades are smoothed out. Plant material originally contains 
the prevailing atmospheric fraction of ^'^C and subsequently, since the material is not 
recycled into the atmosphere, the fraction decreases with the characteristic half life of 

By analysing the ^^C content in the rings of long-lived trees such the Californian 
bristlecone pine, a year-by-year record has been assembled of the cosmic ray flux on 
Earth over the past several thousand years. The data for the last 1000 years are shown 
in Fig. |13| ||^. (In this the following two figures, the vertical axis is oriented so that 
warmer temperatures point upwards.) It is interesting to note that the periods where 
the ^^C deviation exceeds 10 parts per mil correspond to recorded climatic anomalies: 
a) 1100-1250, the so-called Medieval Warm period, b) 1300-1360, the Wolf maximum, 
c) 1460-1550, the Sporer minimum, and d) 1645-1715, the Maunder minimum. The warm 
period that lasted until about 1400 enabled the Vikings to colonise Greenland and wine 
to be made from grapes grown in England. It was followed by a period of about 400 years 
during which the glaciers advanced and cooler, harsher conditions prevailed in most parts 
of the world. 

The Maunder Minimum, when there was an almost complete absence of sunspots, 
would have corresponded to a high cosmic ray flux on Earth and therefore an increased 
cloudiness. This provides a consistent explanation for the exceptionally cold weather of 
this period. Indeed, in every case the excess (c-d) or lack (a-b) of ^^C is consistent with 
the observation of a higher cosmic ray flux leading to more clouds and cooler temperatures, 
and vice versa. 
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Figure 13: History of deviations in the relative atmospheric ^^C concentration from tree- 
ring analyses for the last millennium |jl6| . The data points (dots and open circles) are two 



independent high-precision measurements. The solid line represents a combined fit to a 
large number of other measurements of medium precision. Periods where the deviation 
of the data points exceeds 10 parts per mil are indicated. The first four coincide with 
recorded climatic anomalies. The sharp negative ^^C deviation during the present century 
is the Suess effect, due to the burning of ^^C-depleted fossil fuels. 



These data provide evidence for significant temperature anomalies over the last millen- 
nium that coincided with variations of the cosmic ray intensity. Furthermore they show 
that the Earth has experienced several extended warm and cold spells over the last 1000 
years — with climate swings comparable to the current "anomaly" — whose origin certainly 
cannot be due to anthropogenic greenhouse gases. Whatever caused those earlier spells 
could well be at work today. 



2.3.3 Climatic change since before the last ice age 



The C record has been extended even earlier; the last 10,000 years is shown in Fig. 14 



T7[] and the last 30,000 years in Fig. |T5| [0). There are several important conclusions 



to be drawn from these data. Firstly, there have been very large systematic changes 
in the cosmic ray flux over this period, namely about 40% (where a systematic shift of 
only 1% is known to be associated with a significant climate change, as seen in Fig. 13). 
Secondly, the ending of the last great ice age, 17,000-10,000 years ago, coincided with 
a sharp decrease of the ^"^C fraction (Fig. |1^). Part of this is due to the reduction in 
ice cover and the higher temperatures, thereby increasing the circulation of ^^C-depleted 
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Figure 14: History of deviations of the atmospheric ^'^C concentration over the last ten 
millennia relative to the 1850 value (left-hand scale) [^. The Medieval Warm (MW), 
Wolf maximum (W), Sporer minimum (Sp), Maunder minimum (M) and Suess deviation 
(Su) are indicated. The data points show the Earth's magnetic dipole moment over the 
same period (right-hand scale) ||16||. 
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Figure 15: Variation of the atmospheric ^^C over the last 30 millennia |13]. The two sets 



of lines correspond to theoretical expectations based on changes of the Earth's magnetic 
field; the solid lines assume ±5% of the present-day dipole and the dotted lines, 10%. 



CO2 from the oceans. An estimated decrease in the atmospheric ^^C fraction from these 



effects is 10% [|^. The measured change (Fig. is substantially larger, suggesting that 
this period may have coincided with a decrease of the cosmic flux. 

These large changes in the cosmic ray flux have been caused by changes in the Earth's 
magnetic fieldQ and its associated shielding effect. The geomagnetic field is independently 
measured and found to closely track the cosmic ray flux, as shown in Figs. IT^ and 



15| . This record of changes in the cosmic ray flux, and the expected effects on global 
climate, correlate with the climate records derived from analyses of cores drilled in the 
Greenland ice sheet. However it is important to keep in mind that, even if cosmic rays can 
indeed affect cloud formation, it would also require the existence of suitable atmospheric 
conditions (humidity, aerosol content, etc.) and that these were probably quite different 
during the glacials. 

Finally, on an even longer timescale, it is known that the Earth's magnetic field reverses 
periodically (approximately every million years). During these reversals — which may take 
place fairly rapidly over periods of less than 1000 years — the Earth's field diminishes to 
zero and then re-establishes itself with the opposite polarity. At such times the Earth will 
be subjected to far higher cosmic ray fluxes and the temperatures will be greatly reduced. 
This may perhaps be a factor contributing to the extinction of the dinosaurs 65 million 
years ago, which coincided with a geomagnetic reversal. 

2.4 Conclusions 

In summary, there are two main conclusions to be drawn on the connection between cosmic 
rays and climate change. Firstly, the pattern of systematic change in the global climate 
over the past several thousand years follows the observed changes in the cosmic ray flux 
over the same period; and it is consistent with the explanation that a high cosmic ray flux 
corresponds to more clouds and a cooler climate, and vice versa. Furthermore, the rise of 
about 0.6°C in global temperatures over the last 100 years is consistent in magnitude and 
time dependence with the observed changes in cosmic ray flux — and thereby cloud cover — 
over the same period. A clear and compelling case exists for further study of the effect of 
cosmic rays on cloud formation and climate change — especially in order to improve model 
calculations of the residual effects due to greenhouse gases from the burning of fossil fuels. 

Secondly, since the cosmic ray flux is controlled both by the Sun's magnetic activity 
(solar wind) and by the Earth's magnetic field, it may eventually be possible to make 
long-term (10-1,000 year) predictions of changes in the Earth's climate by detailed mea- 
surements and a deeper understanding of the "magnetic climate" of the Sun and Earth 
and its effect on the cosmic ray flux. In this regard, it will be necessary to understand how 
the Sun changes its magnetic activity and to understand the progression of the Earth's 
magnetic moment. Perhaps even the longer-term Milankovitch cycles are partly influenc- 
ing the global climate through geomagnetic effects; the Earth's magnetic field is due to 
electric currents flowing in the liquid core and in the magnetosphere, and these may have 
been influenced by the periodic gravitational motions of the Milankovitch cycles. 



^The Earth's field is presently weakening (Fig. nj), which may indicate a cooling period ahead. 



3 Cloud formation by charged particles 



Cloud droplets in the troposphere^ form on suitable aerosols above about 100 nm size 
known as cloud condensation nuclei (CCN) |T0|- The new sateUite observations suggest 
that the ions or radicals produced in the atmosphere by cosmic rays may somehow be 
affecting the nucleation, growth or activation of atmospheric aerosols, or affecting the 
creation of ice particles. An increase of CCN concentration could lead to increased droplet 
number concentration, smaller droplet size, less precipitation and therefore longer cloud 
lifetime. If cosmic rays can affect ice nuclei then this would also potentially have a strong 
influence on clouds since the freezing of supercooled water droplets reduces their vapour 
pressure and allows their growth at the expense of neighbouring liquid droplets. 

A microphysical explanation of cloud formation in the presence of ionising radiation 
is presently lacking. In fact there has been very little research on this mechanism and 
present studies do not consider it [|10]. Some experiments with highly charged droplets 
(10~^^-10~^^ C) have, however, shown that they are 10-100 times more efficient at cap- 
turing aerosols than uncharged droplets 



3.1 Droplet formation in a conventional cloud chamber 

Certain aspects of how droplets are induced by ionising radiation are known from the 
development of cloud chambers, but these largely concern high supersaturations of water 
vapour and rapid growth rates. The principle of operation of a cloud chamber can be 



understood from Fig. |T6| which shows the water vapour pressure equilibrium curves for 
small droplets: 
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where p is the vapour pressure, po the saturated vapour pressure at a plane water surface, 
R the gas constant, T the absolute temperature, 7 the surface tension, M the molecular 
weight, p the density, eo the permittivity of free space and r the radius of the droplet. 
The curves divide an upper region of vapour pressure in which water droplets grow by 
condensation (e.g. droplet Di) from a lower region where they shrink by evaporation (e.g. 
D2). If a droplet carries an electric charge Qe then its vapour pressure is reduced by 
increased attraction of the polar water molecules. The growth region is therefore larger 
(Fig. p!6|a) and, in particular, for supersaturations above about 5, all charged droplets 
will grow and reach visible size (e.g. D3). The actual starting point for a droplet is not 
known. However, to set the minimum scale we can estimate the volume occupied by a 
single water molecule as ^7ir^ = ^j/f-, where Na is Avogadro's constant. This indicates 
the effective radius of one water molecule is r ~ 0.2 nm; and a droplet of radius 0.5 nm 
contains about 15 molecules. 



^The troposphere is the lowest level of the atmosphere and the region where there is enough water 
vapour and vertical mixing for clouds to form under suitable conditions. The depth of the troposphere is 
about 6-8 km at the poles, extending to about 17 km over the equatorial regions; it contains about 75% 
of the mass of the atmosphere. There is an overall adiabatic lapse rate of temperature in the troposphere 
by between 6°C (moist air) and 10°C (dry air) per km altitude, reaching a minimum of about -60°C 
at the boundary with the stratosphere (the tropopause). The stratosphere, which extends up to about 
50 km, has a temperature that slowly rises with altitude. This leads to very little turbulence and vertical 
mixing and, in consequence, it contains warm dry air that is largely free of clouds. 
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Figure 16: Formation of charged droplets at 0°C and various supersaturations of water 
vapour, for droplet radii less than a) 5 fim and b) 100 /im. The curves are labelled accord- 
ing to the droplet electronic charge, Qe. The saturated vapour pressure, po, corresponds 
to equilibrium with a plane water surface. Droplets formed in the region above the curves 
will grow (condensation region) and those below will shrink (evaporation region). 
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The necessary supersaturation in a cloud chamber is generated by fast adiabatic ex- 
pansion. For an ideal gas and an adiabatic expansion, 

P^V^ = and 

where 7 is the ratio of specific heats CpjCy (about 1.40 for air and saturated water 
vapour). For example, an expansion ratio V^jVx = 1.3 results in a temperature drop of 
about 30°C and a supersaturation of 4.8. Under these conditions, all droplets with Q > 1 
will grow to diameters of several tens of micrometres in ^ 0.2 s. 

The curves of Fig. |l6|a) indicate how the presence of a small charge facilitates droplet 
formation and growth at high supersaturations, above about 2. The maximum supersat- 
urations in the troposphere are typically less than one per cent above unity which, if this 
were the only mechanism involved, would require large droplet charges for stable growth, 
as indicated in Fig.|16|b). Therefore if ionising radiation is indeed inducing additional 
cloud formation, it must be somehow affecting the aerosols or ice nuclei that efficiently 
seed clouds in the atmosphere. 

3.2 Ions and droplets in the atmosphere 

To provide input for the design requirements of the CLOUD detector, it is useful to esti- 
mate some characteristics of cosmic-ray ionisation and water droplets in the troposphere. 
Some representative parameters at the upper and lower boundaries of the troposphere are 
summarised in Table ^. These regions are well below the ionosphere — which starts above 
the stratosphere at about 80 km altitude — where atoms and molecules are readily ionised 
by solar UV and X-ray photons. 

Table 2: Approximate characteristics of cosmic rays near the top of the troposphere 
(10 km altitude) and at sea level. The maximum and minimum charged cosmic ray 
fluxes, 0, correspond to the sunspot cycle variation at geomagnetic latitudes ^ 40° and 
for energies ^0.1 GeV. The # ion pairs cm~^, riip, are the ionising track densities for 
minimum ionising particles in air. 



Atmospheric characteristics 


Cosmic ray characteristics 


Ah. 


Atmos. 


Pressure 


Density 


Mean. 


Chd. flux, 


Ion pairs 


Recombination 




depth 






temp. 


[m~^sr^-'-s^"^] 


/cm, riip 


coefficient, a 


[km] 


[^int] 


[kN m-2] 


[kg m~^] 


[°C] 


Max. Min. 


[cm^-^] 


^j^3g-lj 


10 


3.0 


26.4 


0.41 


-50 


850 750 


23 


0.64 X 10-12 





11.5 


101.3 


1.22 


15 


100 97 


68 


1.72 X 10-12 



The rate of ion-pair production by cosmic rays, r^p = (pn (100 nip). From Table | 
this rate is about 2x10^ m^^s-^ at sea level and 10^ m-'^s-^ at 10 km altitude. Once 
created, the charged particles will recombine. The electrons are rapidly captured by 
electronegative gases, which become negatively charged as a result. For example, at 0°C 
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and atmospheric pressure, the mean electron attachment time is 190 ns in O2 and 140 ns 
in H2O vapour. 

The ion mobihty is far smaller and consequently so also is the ion recombination rate: 

where n+ and n_ [m~^] are the positive and negative ion densities and a is the recom- 
bination coefficient. Equilibrium is reached when the ion production and recombination 
rates are equal, i.e. when rjp = a n^. Assuming n+ ~ n_ implies the equilibrium ion 

density, n^q = ^Jrip/a. From Table |^, this indicates rieq is about 10^ m~^ at sea level and 
4 X 10^ m~^ at 10 km.Q These estimates apply to still air without aerosols, water droplets, 
etc. The equilibrium ion density rieq oc \/0 so that, for example, the ~12% change in 
cosmic ray flux at 10 km altitude over a sunspot cycle produces ~ 6% change in the ion 
density. 

From equation |l|, the lifetime of an ion is r = l/(a Ueq). This implies ion lifetimes in 
the troposphere of about 6-9 min. In still air, the ions drift vertically in the electric field 
created by the negatively-charged Earth and the positively-charged ionosphere. The field 
strength is E' ~ 100 V/m at sea level, producing an ion drift velocity of about 0.022 m s~^. 
This results in an ion drift distance of about 12 m before recombination. 

In the CLOUD detector, in the absence of a drift field, the ions will diffuse a projected 
distance ax = y/2Dt in a time t, where D is the diffusion coefficient. For water vapour, 
D = 1.6 X 10~® m^s~^ at conditions corresponding to sea level and 4.1 x 10~^ m^s~^ at 
those corresponding to 10 km altitude, indicating ion diffusion distances of 4-6 cm before 
recombination. 

In order to be sensitive to long growth times in CLOUD, the droplets must not be 
removed by gravitational fall before observation. This sets a maximum droplet size for 
studying long growth times, which can be estimated as follows. Stoke's formula for the 
viscous force F on a sphere of radius r moving at a velocity v is 

F = —6'K'qrv, (2) 

where rj is the dynamic viscosity (1.7 x 10^^ Pa s for air at 0°C). This indicates a terminal 
velocity Vo = {2r'^ pg) / {9ri) , or Vo [/um s^^] = 130 [/im^]. A droplet of radius 1 pm 
therefore has a terminal velocity 130 /xm s~^ and falls by about 8 cm during a time of 
10 min. 

Very small droplets will be suspended by Brownian motion. We can estimate their 
maximum size as follows. The projected Brownian displacement is = y/2Dt, where 
D = iJ,kT, and k is the Boltzmann constant. The mobility is = Vdrift/F = l/{67Trir) 
(from equation where F is the force causing the drift. Substituting gives the rms pro- 
jected displacement of a droplet due to Brownian motion, ax = \J {kTt) / (Snrir). Brownian 
motion will dominate over gravitational fall when the projected Brownian displacement 
is about equal to the droplet radius r in the time the droplet would have taken to fall a 
distance r. This indicates that Brownian displacement is comparable to gravitational dis- 
placement at a droplet radius of about 0.7 /xm, and below this value droplets will remain 
suspended for extended periods. 

^For comparison, the droplet density in clouds is 10^-10® m"'^ over oceans and 10^- 
land. 
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In summary, in order to allow for the possibility of long formation times in the atmo- 
sphere, the CLOUD detector should have the capability to observe droplet growth over 
periods of up to about 10 min, i.e. comparable to the ion lifetime, and the dimensions of 
the chamber should allow for a fiducial volume that is approximately 10 cm away from any 
boundary. This indicates a chamber dimension of about 50 x 50 x 50 cm^ is appropriate. 



4 The CLOUD experiment 

The basic purpose of the CLOUD detector (Fig. |T^) is to confirm, or otherwise, a di- 
rect link between cosmic rays and cloud formation by measuring droplet formation in a 
controlled test-beam environment. The CLOUD detector is essentially a cloud chamber 
that is designed to duplicate the gases, water vapour, aerosols and temperature/pressure 
conditions found in the Earth's troposphere. 

Since the microphysical mechanisms responsible for cloud formation by cosmic rays 
are poorly understood, it is necessary for the experimental tests to cover a broad range 
of variables and conditions. For instance, it is not yet even known at which altitude 
(i.e. temperature and pressure conditions) the clouds are being formed. Also, since 
the maximum water vapour supersaturation in the troposphere is only about one per 
cent, droplet growth times may be quite long, as discussed in the previous section. The 
combined effects of ionisation and aerosols/ice particles is probably important (Section^. 
The ionisation of existing water droplets by charged particles requires study; charged 
droplets have a lower vapour pressure and will therefore grow at the expense of uncharged 
droplets. It is important to measure the sensitivity to a range of track ionisation densities, 
dE/dx, and ion pair concentrations that duplicates those of cosmic rays. This indicates 
the need for beam particles ranging from protons to heavy ions. Finally, an investigation 
should be made of the possibility that droplets are formed in relatively "rare events" , such 
as very highly charged ions or droplets, perhaps due to heavy nuclei. 

Cloud chamber data under these conditions have never been previously obtained. 
Although C.T.R. Wilson's cloud chambeiQ was extensively developed in the first half 
of this century, it was optimised for experimental particle physics and operated under 
conditions far removed from those of the troposphere. 



4.1 Detector 



The CLOUD detector (Fig. comprises the following main components: 



A temperature-controlled cloud chamber that can be operated in the temperature 
range -60°C to 30°C and at pressures from a vacuum to 1.5 atm. The upper pressure 
limit is to allow for adiabatic expansions down to a final pressure of 1 atm. An 
electrode structure (field cage) in the cloud chamber provides a simple clearing field 
for measurements with low residual ionisation. 



''it is interesting to note that C.T.R. Wilson had the inspiration for the cloud chamber while observing 
meteorological phenomena on the mountain of Ben Nevis in 1894. (The phenomena were not particle 
tracks, however, but "coronas" around the Sun and glories, where the Sun glows around shadows in the 
mist.) 
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Figure 17: Schematic diagram of the CLOUD chamber. 
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• A surrounding liquid bath enclosed by super-insulation to provide stable and precise 
temperature conditions for the cloud chamber. A temperature stability of about 
0.01 K is required since this is equivalent to a change of water vapour supersaturation 
by about 0.05%. The investigation of long growth times for the droplets requires 
that the operating conditions be held steady for long periods. The temperature is 
adjusted by heat-exchange coils immersed in a suitable liquid coolant in the bath and 
connected to a He compressor. The temperature of the circulating He is controlled 
by temperature monitors immersed in the coolant. 

• A piston to provide adiabatic expansions for investigating rapid droplet growth 
times. Small piston expansions and compressions also provide a technique for fine 
control of the water vapour supersaturation in the cloud chamber. 

• An optical system comprising illumination and stereo microscope cameras (recu- 
perated from old bubble chamber equipment) which can be operated with film or 
CCD readout. The latter is to allow direct digitization of the data for analysis. 
The illumination and optical system must be capable of detecting and measuring 
water droplets of sizes down to a minimum of about 0.5 fim, i.e. the smallest that 
can scatter light and therefore constitute a cloud. The beam pipe connecting to 
the cloud chamber is equipped with further detectors to measure the ion/aerosol 
concentration, size distribution and charge. 

• Temperature and pressure monitors. 

• A gas system providing air with adjustable amounts of water vapour and aerosols. 
A mass spectrometer is required for physical and chemical analysis of the incoming 
and outgoing gas/aerosol mixture. 

• A scintillation counter hodoscope to measure the incident beam. The beam must 
be spread over a large area to simulate the quasi-uniform cosmic radiation. 

• A DAQ system. 

4.2 Experimental goals 

The goals of the proposed study are threefold: 

1. Meaisurements of water droplet (cloud) formation by ionising radiation, 

according to the following variables: 

• Temperature and pressure conditions. 

• Supersaturation of water vapour. 

• Presence of aerosols and trace condensable vapours. 

• Presence of frozen droplets. 

• Presence of existing water droplets. 

• Particle flux and ionisation density (p, heavy ions. . . ). 

• Time for droplet growth. 
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2. Understanding of the microphysical processes involved in cloud formation by 
ionising radiation. This wiU involve Monte Carlo simulation of the processes and a 
comparison with the CLOUD results. 

3. Estimation of the global cloud cover due to cosmic rays, using atmospheric 
data recorded by satellite experiments. 

5 Conclusions 

The observation of an apparent connection between cosmic rays and global climate offers 
a unique opportunity for particle physics to make a major contribution to the problem 
of global warming — at a relatively modest cost. The proposed CLOUD experiment will 
unambiguously confirm, or otherwise, a direct link between cosmic radiation and cloud 
formation, and allow a quantitative understanding of the physical processes involved. 
This, in turn, would provide an unambiguous and physically explicable connection be- 
tween solar and climatic variability — representing a great advance in the understanding 
of the relationship between the Earth and its star. 

Although this experiment primarily addresses atmospheric science, it is clearly in the 
domain of particle physics and of CERN's expertise and experimental facilities. The 
CLOUD experiment will have the immediate impact of potentially resolving one of the 
important unknown effects that may have so far prevented reliable calculations of global 
warming from greenhouse gases — an issue of profound economic and social importance to 
the world. Furthermore, if the link is confirmed, then it may eventually become possible 
to make long-term predictions of systematic changes in the Earth's climate, by monitoring 
and predicting the evolution of the cosmic ray flux and the "geomagnetic climate" of the 
Sun and Earth that modulates its intensity. 
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